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1 Introduction 



With the evidence for neutrino masses and mixing in solar pQ and atmospheric |2J neutrino 
experiments, studies of the lepton sector have gained importance as a path to physics 
beyond the Standard Model. The most elegant and widely accepted explanation for small 
neutrino masses is provided by the see-saw mechanism [3], in which a large Majorana mass 
scale Mr of right-handed neutrinos drives the light neutrino masses down to or below the 
sub-eV scale, as required by the experimental evidence. A priori, the fundamental scale Mr 
can be of the order of the GUT scale, and may thus be unaccessible for any kind of direct 
experimental tests. However, neutrino mixing implies lepton-flavor violation (LFV), which 
is absent in the Standard Model and provides indirect probes of Mr. While lepton-flavor 
violating processes are suppressed due to the small neutrino masses if only right-handed 
neutrinos are added to the Standard Model jlj, in supersymmetric models new sources of 
LFV exist. For example, virtual effects of the massive neutrinos affect the renormalization 
group equations (RGE) of the slepton mass and the trilinear coupling matrices, and give 
rise to non-diagonal terms inducing LFV. 

Assuming the experimentally favored large mixing angle (LMA) MSW solution of the solar 
neutrino anomaly, one can expect lepton-flavor violating /i and r decays with branching 
ratios close to the current experimental bounds [3]. Some of the existing bounds will be 
improved significantly in the near future. The current experimental limits (future sensi- 
tivities) on low-energy lepton-flavor violating processes involving charged leptons can be 
summarized as follows: 

• Br(fi -> erf) < 1.2 • Kr n (l(r 14 ) Ej 

• Br{r -> erf) < 2.7- KT 6 [8 

• Br{r -> yq) < 1.1 ■ Kr 6 (l(T 9 ) |E| (1) 

• Br(n + -> e + e + e~) < 1.0 • 1(T 12 [TT] 

• Rifi'Ti -> e~Ti) < 6.1 • 10" 13 (10" 14 ) [HI US! 

Here, the observable R denotes the cross-section normalized to the total muon capture rate. 
The MECO experiment aims at a sensitivity for pTAl — > e~Al below R « 1CT 16 14]. In 
the farther future, the PRISM project plans to provide beams of low-energy muons with an 
intensity increased by several orders of magnitude, so that it may become possible to reach 
Br (ji -> erf) « KT 15 [15 , Br -> e+e+e") « KT 16 [16] and R(fU~Ti -> e~Ti) ta 10" 18 
[T7j (see also the review |18j). Searches for r — >• /i7 at the LHC or SUPERKEKB are 
expected to probe LFV in this channel at the level of Br w 10 -9 [TUj . 
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The above processes in the context of supersymmetric see-saw models have been considered 
in several previous studies (see e.g. 113120112111221 I2H 121] )• In [H] it has been pointed out 
that the corresponding branching ratios and cross-sections exhibit a quadratic dependence 
on the right-handed Majorana neutrino mass scale Mr. Therefore, the exploration of these 
processes provides very interesting possibilities to constrain Mr. In the present paper, we 
sharpen the current knowledge of these constraints by investigating in more detail which 
information about the right-handed Majorana masses can be extracted from measurements 
of the processes (JIJ. It is assumed that the right-handed Majorana masses are degenerate at 
the scale Mr. We focus on the recently proposed post-LEP mSUGRA benchmark scenarios 
[25] . and take into account the uncertainties in the neutrino parameters. In addition, we 
show by how much the sensitivity to Mr will improve with future more precise neutrino 
data. Our work updates and extends previous studies in several directions. Firstly, the 
mSUGRA scenarios of j2H] have been developed particularly for linear collider studies, but 
have not yet been applied to lepton-flavor violating processes at low energies. Our study 
clarifies the model-dependence of the latter for this very relevant set of mSUGRA models. 
Secondly, the neutrino input in our analysis is varied in the ranges allowed by present data. 
The results are compared to expectations for more precise neutrino measurements in the 
future. Thirdly, we consider degenerate as well as hierarchical neutrino spectra, and study 
the impact, a future determination of the absolute neutrino mass scale would make. Finally, 
following [2B] it is demonstrated that the influence of the mSUGRA scenarios in the tests of 
Mr can be reduced by normalizing Br{li — > Ijj) to the corresponding SUSY contribution 
to the muon anomalous magnetic moment. 

The paper is organized as follows. In section 2, we discuss the supersymmetric see-saw 
mechanism and the renormalization group evolution of the neutrino and slepton mass ma- 
trices. In section 3, the rare decays Zj — > Ijj, \x — > 3e as well as /i-e conversion in nuclei 
are briefly reviewed, and the most important results for our investigations are displayed. 
Also the anomalous magnetic moment of the muon and the correlation with U — > Ij'j is 
discussed there. Section 4 summarizes the input parameters of the mSUGRA benchmark 
scenarios and the experimental neutrino data used in the analysis. The numerical results 
of our studies are presented in section 5, and conclusions are drawn in section 6. 
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2 Super symmetric see-saw mechanism 

The supersymmetric see-saw mechanism is described by the term 

1 



W„ = --v%Mv c R + v%Y v L-H 2 



(2) 



in the superpotential, where VR a [a = e, ji, r) are the right-handed neutrino singlet fields, 
L a denote the left-handed lepton doublets and H 2 is the Higgs doublet with hypercharge 
+ 2- The 3x3 matrix M is the Majorana mass matrix, while Y v is the matrix of neutrino 
Yukawa couplings leading to the Dirac mass matrix mc = Y V {H 2 ) , (H 2 ) = f sin/5 being 

(H°) 

the H 2 vacuum expectation value with v = 174 GeV and tan/3 = |gSy- Light neutrinos can 
be naturally explained if one assumes that the Majorana scale Mr of the mass matrix M 
is much larger than the scale of the Dirac mass matrix mo, which is of the order of the 
electroweak scale. At energies much smaller than Mr one has an effective superpotential 
with 

(3) 



W e Jf = -{Y U L ■ H 2 ) T M-\Y U L ■ H 2 ). 

The corresponding mass term for the left-handed neutrinos v^ a is then given by 

1 



^\M v v L + h.c, 



(4) 



where the mass matrix 



M v = mlM^rno = Y^M~ x Y v {v sin (3f 



(5) 



is suppressed by the large Majorana scale Mr. In the following we work in the basis 
where the charged lepton Yukawa coupling matrix YJ 1 and the Majorana mass matrix M 
of the right-handed neutrinos are diagonal, which is always possible. The matrix M v is 
diagonalized by the unitary MNS matrix U, 



U T M U U = diag(mi, m 2 , m 3 ), 
that relates the neutrino flavor and mass eigenstates: 



(6) 



u 



v 2 



(7) 



1 Therefore, we do not have to discriminate flavor and mass eigenstates for charged leptons, i.e. l e ,fj,,r = 
'1,2,3 = e,fi,T 
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In general, U can be written in the form 

U = \/-diag(e^ 1 ,e^ 2 ,l), (8) 
where <fti, <p2 are Majorana phases and V can be parametrized in the standard CKM form: 
/ C13C12 C13S12 s 13 e~ i<p \ 



V 



-C23S12 - ^235130126^ C23C12 - S 2 3Si 3 S 12 e ltp S 23 Cl 3 



(9) 



\ S23S12 — C23Si3Ci2e* ¥ ' — S23C12 — C23Si3Si2e* v C23C13 J 

The experimental data on neutrino oscillations determine or at least constrain the mixing 
matrix V and the differences of the squared mass eigenvalues m ; at a scale not far from 
the electroweak scale. We will therefore identify these two scales in our analysis. Using 
the results of recent neutrino fits and making some further necessary assumptions on the 
neutrino spectrum one can reconstruct M v {Mz) from (JHJ). 

2.1 Renormalization group evolution of the neutrino sector 

In order to calculate the lepton-flavor violating contributions to the slepton mass matrix in 
a top-down approach from the unification scale Mx ~ 2 ■ 10 16 GeV to the electroweak scale, 
we first need to evolve the neutrino mass matrix M u (Mz) to Mx- Below Mr, the one-loop 
RGE in the MSSM is given by [T7j 

±M V = ( (-6g 2 2 - %g\ + Tr^Yu)] M v + ((Y?Yi)M v + M^Y^f)) (10) 



dt v 167T 2 VV ^ 5^ 

with the U(l) and SU(2) gauge couplings g\ and g 2 , and the Yukawa coupling matrices 
Y\j and Y\ for the charge |-quarks and charged leptons, respectively. The corresponding 
evolution equations for g\^ Yy and Yi can be found in [2E|- The RGE is linear in M v and 
can thus be solved analytically |2*7j : 

M v (t)=I(t)-M v (0)'I(t), t = ]n(jj~y (11) 

Since the evolution is dominated by the gauge and third generation Yukawa couplings one 
obtains, to a good approximation: 



with 



I(t) = / s / t diag(l,l,/ T ) (12) 



w = ex p(i^i <( " 3fe2 "^ i2) ^) (13) 

I t {t) = expf J- f Z\Y t \ 2 dA (14) 



167T 2 JO 



I T (t) = exp finiW) . (15) 
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To calculate these factors, the MSSM RGEs for the gauge and Yukawa couplings are solved 
in one-loop approximation, neglecting threshold effects. 

In order to proceed with the evolution from Mr to Mx we use directly the matrix Y u of the 
neutrino Yukawa couplings. From (j3J) and (jUJ) one finds (201 

Y u = — L-diag(*/Mi, \[m 2 , Jm~ 3 ) ■ R • diag(v^r, ^m~ 2 , 0^) • U\ (16) 

where Mj are the Majorana masses of the right-handed neutrinos and R is an unknown 
orthogonal matrix. As we will see, the lepton-flavor violating terms in the slepton mass 
matrix depend on Y v only through the combination YJY V . In this work, we assume the 
right-handed Majorana masses to be degenerate at Mr (Mi = M 2 = M3 = Mr) and the 
matrix R to be real. Then the product YJY V simplifies to 

Mr 

Y l Y v = o • 2 a U ' diag(mi, m 2 , m 3 ) • W, (17) 
v z sin p 

thus being independent of R. Therefore this class of models is highly predictive and often 
used for phenomenological studies. In addition, one obtains more conservative upper bounds 
on Mr for real R because a complex matrix R generically leads to larger values of Y^Y V 
and thus to larger branching ratios Br(li — > Ijj) as shown in [20]. Furthermore, since the 
Majorana phases 4>i an d 02 defined in (JSJ) also drop out in (|I7j). U can be replaced by V in 
f|16|) . Finally, the neutrino masses rrii and are evaluated from M v at Mr using (jUJ). The 
resulting matrix Y u (Mr) is then evolved from Mr to Mx using the one- loop RGE (20] 

and keeping the product YJY V on the r.h.s. of ()18j) fixed at Mr. The running of the right- 
handed mass matrix M between Mr and Mx is negligible, as we have checked numerically. 

2.2 Renormalization group evolution of the slepton sector 

Having evolved the neutrino Yukawa couplings or, more specifically, the product YjY v , 
to the unification scale Mx, one can now run the slepton mass matrix from Mx to the 
electroweak scale assuming the mSUGRA universality conditions at Mx' 

ml = mil, m 2 R = m 2 l, A = A Y h (19) 

where mo is the common scalar mass and Aq is the common trilinear coupling. For the 
present analysis, we adopt the mSUGRA benchmark scenarios proposed recently in [23] for 
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linear collider studies. The charged slepton (mass) 2 matrix has the form: 

(771- ( 771- V \ 
It ^ l,n> \ (20) 
mf mf I 
Ilr Ir ' 

where m~ , m~ and m~ are 3x3 matrices, m~ and m 2 being hermitian. The matrix 

lh l R l LR t-L lR 

elements are given by 

{m} L )ab = (m 2 L ) ab + 5 ab (m 2 a + m|cos(2/3) (-i + sin 2 ^)) (21) 
(mfjab = (m 2 R ) ab + 6 ab (m 2 a -m|cos(2/5) sin 2 9 W ) (22) 
( m I LR )ab = A ab vcosp - 5 ab mi a fxtanp, (23) 

9w being the Weinberg angle and \x the SUSY Higgs-mixing parameter. After evolution 
from Mx to Mz, one has 

m\ = mil + (Sm 2 L ) USSM + 5m 2 L (24) 

m 2 R = mil + (5m 2 R ) MSSM + Sm 2 R (25) 

A = A Q Yi + 5v4 MSSM + 5 A, (26) 

where (Sm 2 ^ R ) MSSM and (5 A) MSSM denote the usual MSSM renormalizat ion-group corrections 
j2Hl which are flavor-diagonal. In addition, the presence of right-handed neutrinos radia- 
tively induces flavor off-diagonal terms denoted by 8m 2 LR and 5 A in (J24j) to (J26j) . These 
corrections are taken into account in the approximation ^Hj 

Sm 2 L = -_L(3mg + X2)(^n)ln(^) (27) 
Sm 2 R = (28) 

It is these terms which give rise to lepton-flavor violating processes such as U —>■ Iff and 
/i-e conversion. 

The physical charged slepton masses are then found by diagonalizing (J2(Jj) using the 6x6 
unitary matrix Up 

Utm-Uf = diag(m 2 i , ...,m~., (30) 

Correspondingly, the slepton mass eigenstates are expressed in terms of the gauge eigen- 
states by 

U = (Uf) ai l La + (Uf)( a +3)iliia, i = 1, — ,6; a = e,fi,r. (31) 
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Similarly to (J21|) . the 3x3 (mass) 2 matrix of the SUSY partners of the left-handed neutrinos 
is given by 

(ml) ab = (m 2 L ) ab + -S ab m 2 z cos(2/3), (32) 

where m 2 L can be taken from (J24|) . The partners of the right-handed neutrinos are very heavy 
and can therefore be disregarded. After diagonalization with the unitary 3x3 matrix U$, 

Ulm%U = diag^-^m^m-J, (33) 

the mass eigenstates z/j are related to the gauge eigenstates by 



/ 




\ 
















V 




J 







3 LFV low-energy processes and — 2 

3.1 The radiative decays U — > ^7 

The effective Lagrangian for l~ — > /~7 is given by 

£eff = faapF* (AlP L + A%P R ) k, (35) 

where F a/3 is the electromagnetic field strength tensor, = ~ [j a , jp] and Pr l = |(1±7 5 ) 
are the helicity projection operators. The coefficients A l [ R are determined by the photon 
penguin diagrams shown in Fig. 1 with charginos/sneutrinos or neutralinos/charged sleptons 
in the loop. 
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From (|35|) one obtains the following decay rate for l { — > £ • 7 [5]: 

r (C - /Jt) = f ^ (l^, + ^l 2 + \A% + Al\ 2 ) . (36) 

The superscript c (n) refers to the chargino (neutralino) diagram of Fig. 1, while the flavor 
indices are omitted. Because 3> and m 2 ^ is diagonal (see (J23)) and (|2EJ)), one 
has A/j A L [TH], [2D]- The dominant amplitudes in (}3T)j) are approximately given by 



327T 2 y^mvi/ cos /3 a=1 



- fe=i v fc 

1 / 
x 



(-3 + 4r^-(r= fc ) 2 -21nr^) (37) 



^ ^ -i^ 2 tan^^5:5-(C' i v) ol ((C> 7 v) a2 + (C> i v) al tan^ 



32tH a=l^i< 



xiU-OAUfh+W-, -3 1 " KkV + 2C fc lnC fe (38) 



with 



(39) 



the chargino diagonalization matrices Ol, Or and the neutralino diagonalization matrix 
On. The mass eigenvalues of the charginos and neutralinos are denoted by m^- and m^o, 
respectively. The numerical calculations discussed later are performed with the full expres- 
sions for A c £ n and A ^ 1 , which can be found in |3] and [3P] . 

Note that there is no difference between the rates of — > Ij^y and if — > if 7 at the one-loop 
level and no CP violating observables can be constructed at this level of perturbation theory 
[301 EI]. We therefore do not distinguish between Br(l^ — > Ijj) and Br (If — > ifj) in the 
following. 



3.2 Br (/i 3e) and R (fi~N -> e~iV) 

The processes /i — > 3e and /i~iV — > e~iV are dominated by photon penguin contributions. 
As a consequence, one has the following model-independent relations 19 : 

BrQ.~3e) „ £ 8 U . U\ „ ? , ^ (4Q) 
5r(/x -> e7) 8vr 3 \ m 2 4 / v y 

« 6 ■ 1(T 3 for Titanium, (42) 
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where is the total decay width of the muon, F(q 2 ) is the nuclear form factor and Z (Z e ff) 
is the electric (effective) charge of the nucleus. From the above and (0) one can see that 
the present experimental upper limits on Br (fi —>■ 3e) and R {fi~N — > e~N) constrain LFV 
considerably less than the current limit on Br (n — > ej). However, a future measurement 
of R in the range of R ([i~Ti — > e~Ti) « 10 -18 ^7] as mentioned in the introduction would 
provide a more sensitive test than the corresponding future sensitivity Br (/i — > e^y) ~ 10~ 15 . 
We have examined the relations (}4*0"|) to (}4*2"|) numerically for the neutrino parameters and 
mSUGRA scenarios presented in the next section using complete analytic expressions [5], 
and have found 



Br(fi — > 3e) 
Br(n — > ej) 
R(fi~Ti —>■ e~Ti) 
Br(fi — > ej) 

in good agreement with the above estimates. 



(6 - 7) ■ 1(T 3 (43) 
(5 - 7) • 1(T 3 (44) 



3.3 Anomalous magnetic moment of the muon 

The supersymmetric contribution to \{g^ — 2) is described by the diagrams of Fig. 1 for 
i = j = 2. In this case the effective Lagrangian (j33j) yields 

K = ^ K + Af) (45) 

and, with (f3lj|). the relation j2H] 

Br (lj -> fj-7) ^ a mf. 

I^l 2 T im 2 

where Tj denotes the total decay width of lepton ^ and A 1 ^ ^> A l [ has been used. It will be 
shown later that this ratio varies less with the SUSY parameters and thus provides a less 
model-dependent test of LFV than Br (Zj — > Ijj) alone. 



A 13 



Af + Af 



(46) 



4 Input parameters 

4.1 mSUGRA benchmark scenarios 

In this paper we focus on the mSUGRA benchmark scenarios proposed in |25j . The theo- 
retical framework of these scenarios is the constrained MSSM with universal soft supersym- 
metry breaking masses and imparity conservation. Sparticle spectra corresponding to these 
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scenarios are consistent with all experimental and cosmological constraints, in particular 
with 

• direct sparticle searches 

• b — > 57 

• cosmological relic density, with the lightest neutralino as lightest SUSY particle and 
dark matter candidate 

• Higgs searches 

This class of models involves five free parameters: the universal gaugino mass m.1/2 and the 
universal scalar mass tuq at the GUT scale, the ratio tan /3 of the Higgs vacuum expectation 
values, the sign of the Higgs mixing parameter /x, and the universal trilinear coupling 
parameter A . The values of these parameters for the benchmark scenarios are listed in 
Tab. A is set to zero in all scenarios. Further details can be found in [25 . 



Scenario 


m 1/2 /GeV 


m /GeV 


tan/? 


sgn(/i) 




A 


600 


140 


5 


+ 


2.8 


B 


250 


100 


10 


+ 


28 


C 


400 


90 


10 


+ 


13 


D 


525 


125 


10 




-7.4 


E 


300 


1500 


10 


+ 


1.7 


F 


1000 


3450 


10 


+ 


0.29 


G 


375 


120 


20 


+ 


27 


H 


1500 


419 


20 


+ 


1.7 


I 


350 


180 


35 


+ 


45 


J 


750 


300 


35 


+ 


11 


K 


1150 


1000 


35 




-3.3 


L 


450 


350 


50 


+ 


31 


M 


1900 


1500 


50 


+ 


2.1 



Table 1: Input parameters of the mSUGRA benchmark scenarios and the predicted shift 
Sa^ in |(^-2) |2S1. 

Also given in Tab. ^is the corresponding shift in the muon anomalous magnetic moment. 
The current 1.6 o discrepancy between the measurement of — 2) and the Standard 
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Model prediction [32] amounts to 

5a^ = (25 ± 16) • 1(T 10 . (47) 

Scenarios with relatively light sparticle masses below 500 GeV (e.g. B, C, G, L) are in 
better agreement with the above value of Sa^ than scenarios with heavier sparticles (e.g. E, 
F, H, M). Moreover, (jTfj) favors a positive sign for //. 

4.2 Neutrino data 

Solar and atmospheric neutrino experiments provide clear evidence for neutrino oscillations. 
The favored interpretation of the experimental results on solar neutrinos suggests v e —>■ 
oscillations driven by the mass squared difference AmJ 2 = m\—m\ in the range of the LMA 
solution, while the results on atmospheric neutrinos are interpreted by — > v T oscillations 
driven by Am^ = m\ — m\ in the case of three active neutrinos. For the present analysis, 
we use the global fits in a three-neutrino framework performed in (33] • in [HI] it has been 
pointed out that the inclusion of the SNO result in a two-neutrino analysis of the solar 
neutrinos implies only minor changes. 

We also consider the improvement in our knowledge of the neutrino sector expected from 
future neutrino experiments, and discuss the consequences of these future accomplishments 
for the tests of the Majorana scale Mr considered in this paper. We always assume the 
present best fit values of the neutrino parameters to remain unchanged. The future im- 
provements in the experimental errors of these parameters anticipated for a perspective 
view are summarized below together with other relevant expectations: 

• Am\ 2 an d s i n2 2#i2 ■' The long-baseline reactor experiment KAMLAND is designed to 
test the LMA MSW solution of the solar neutrino problem. Data taking is expected to 
start in 2002 and the solar neutrino parameters will be determined with an accuracy 
of 5(Am 2 12 )/Am 2 12 = 10% and <5(sin 2 2# 12 ) = ±0.1 within three years of measurement 

m 

• Am 2 3 and sin 2 2# 23 ■' The atmospheric oscillation parameters will be determined by the 
long-baseline accelerator experiment MINOS with an accuracy of <5(Am 2 3 ) / 'Am| 3 = 
30% and 5(sin 2 2fl 23 ) = ±0.1 [H]. 

• sin 2 2# 13 : The CHOOZ reactor experiment restricts the angle #13 to sin 2 2# 13 < 0.1 
|37j . The long baseline experiment MINOS [HE] can probe the range sin 2 26*13 > 0.02- 
0.05. A future superbeam, a neutrino factory [3H] or the analysis of the neutrino 
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energy spectra of a future galactic supernova may provide a sensitivity at the 
level of a few times 10~ 3 to 10~ 4 . To explore the potential of future neutrino studies 
we take 5(sin 2 20 13 ) = 3 ■ 1(T 3 . 

• The neutrino mass spectra: The inverse hierarchical spectrum with two heavy and 
a single light state is disfavored according to a recent analysis [40 j of the neutrino 
spectrum from supernova SN1987A, unless the mixing angle #13 is large (compare, 
however, |Hj). We therefore restrict our analysis to the direct (normal) hierarchy. 
LFV rates for inverse hierarchical schemes lie in the intermediate range between the 
extreme cases we discuss, Br (degenerate) -C Br (inverse) < Br (hierarchy), as pointed 
out in j2U- 

• The Dirac CP phase (p: Even at a neutrino factory, one will only be able to distinguish 
ip = from it/2 if Am 2 2 > 10~ 5 eV 2 . For this reason we vary (p in the full range 
< <p < 2vr E2|. 

• The neutrino mass scale: While neutrino oscillation experiments provide information 
on the neutrino mass squared differences Amfj, the absolute scale of the neutrino 
masses is not known so far. Upper bounds can be obtained from the neutrino hot dark 
matter contribution to the cosmological large scale structure evolution and the Cosmic 
Microwave Background, from the interpretation of the extreme energy cosmic rays in 
the Z-burst model, tritium beta decay experiments and neutrinoless double beta decay 
experiments [12] • A next generation double beta decay experiment like GENIUSI, 
MAJORANA, EXO, XMASS or MOON will test the quantity m ee = | Vge* 2 *^] 
down to 1CT 2 eV. Since = sin 2 26 l i 3 /4 < 0.025, the contribution of m 3 drops out 
and a bound m ee < 10 -2 eV will imply mi < 10~ 2 eV/ cos 2#i2- If one further assumes 
that KAMLAND measures sin 26 12 with 5(sin 2 2# 12 ) = ±0.1, one obtains the bound 
mi < 3 • 10 -2 eV. On the other hand, a large mass mi could be tested by future 
tritium beta decay projects. A positive signal at the final sensitivity of the KATRIN 
experiment would imply mi = 0.3 ± 0.1 eV jH]. Such a value would be compatible 
with the recent evidence claim for neutrinoless double beta decay |4*o] . 

For the present purposes, typical hierarchical and degenerate neutrino mass spectra are 
parametrized as follows: 

(a) hierarchical vi and degenerate Ur\ 








(48) 



13 



.\/| M, M 3 Mj 



(49) 



(b) quasi-degenerate ul and degenerate Vr 



mi, 
M 1 -■ 



m 2 
M 2 



mi + 
M 3 = 



1 A 2 1 2 
— Am 12 , f»3 »s mi H Am 23 



2m! 

Mr 



2mi 



(50) 
(51) 



where mi ^> i/Ami ^> 



Am 2 



12- 



The product of Yukawa couplings YJY V appearing in the renormalization group corrections 
to the left-handed slepton mass matrix (|2*7)l can then be approximated by 



(b) (YX 



ab 



ab 



v 2 sin 2 (5 

M R 

v 2 sin 2 j3 



Am 



\ Am 



?K 2 n*2 + KaKs 

23 



mi$ab + 



Am 23 ( Ami 
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2m 1 VAmlg 2 



a3 Vfo3 



In both cases the largest branching ratio for — > L7 is expected in the channel r 



because of | V33 | 



V 32 V£\ and Am 2 q > A 



2:1 



ni 



12- 



The decays /x — > and r 



(52) 

(53) 

■> ^7 
-> e7 



are suppressed by the smallness of Am 2 2 and V13. In the case (b), there is an additional 



suppression by yAm 23 /mi or y Am 2 2 /mi relative to the case (a). 

In the following analysis, the neutrino parameters are varied in the ranges specified in 
Tab. [21 characterizing the present knowledge and future prospects. 



parameter 


best fit value 


present 


future 


tan 2 6*23 


1.40 


+1.64 
-1.01 


+1.37 
-0.66 


tan 2 #i3 


0.005 


+0.050 
-0.005 


+0.001 

-0.005 


tan 2 6*12 


0.36 


+0.65 
-0.16 


+0.35 
-0.16 


Am 2 2 /10" 5 eV 2 


3.30 


+66.7 
-2.3 


+0.3 
-0.3 


Am^/10" 3 eV 2 


3.10 


+3.0 
-1.7 


+1.0 
-1.0 




to 2vr 


mi/eV 


to 0.03 (0.3ljnJ) 



Table 2: 90% CL fits of neutrino parameters characterizing the present and future un- 
certainties. The range of the neutrino mass scale mi refers to a hierarchical (degenerate) 
spectrum. 
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Note that all parameters are varied simultaneously. The Majorana mass scale Mr is treated 
as a free parameter. This contrasts with other approaches ^HJ 120] where Yukawa coupling 
unification 

|Y„ 3 |=Y t atM x (54) 

is assumed, \Y V3 \ 2 being the largest eigenvalue of YjY u . Fig.|2]shows the normalized Yukawa 
coupling |y, 3 |/y at Mx as a function of the Majorana mass Mr. One can see that the 
assumption would fix the Majorana mass scale to Mr « 4 • 10 14 GeV for hierarchical 
neutrinos and to Mr « 7 • 10 13 GeV in the degenerate case. 

Fig. |21 also shows that for large values of Mr, the Yukawa coupling \Y V2 \ eventually gets too 

\y i 2 

strong for perturbation theory to be valid. Therefore, we restrict \Y V3 \ to values ^ < 0.3, 
which implies the consistency limits Mr < 2 ■ 10 15 GeV in the hierarchical and Mr < 
3 • 10 14 GeV in the degenerate case. It should be stressed in this context, that since the 
negative mass shift 5m\ given in (|27|) is driven by the neutrino Yukawa couplings, the 
slepton masses decrease with increasing Mr. We have checked that in the perturbative 
region of Y Vs defined above, the slepton masses do not violate existing lower mass bounds, 
in particular the LEP bound > 81 GeV [H|. 

5 Numerical Results 

The dependence of Br(li — > Ijj) on the right-handed Majorana mass scale Mr is displayed 
in Figs. 3, 4 and 5 for those two mSUGRA scenarios listed in Tab. ^ which lead to the 
largest and smallest branching ratios. The sensitivity on Mr for all benchmark mSUGRA 
scenarios defined in Tab. Q is summarized in Tab. 3. The present bounds on Br{r — ► ej) 
and Br(r — > fi'j) set relatively weak constraints on Mr and are therefore not included in 
Tab. 3. For each scenario, the neutrino input is varied in the ranges allowed by present 
and/or future experiments. 

Comparing Figs. 3 and 4 with Fig. 5, one can see that Br(fi — > ej) is more strongly 
affected by the uncertainties in the neutrino parameters than Br(r — > fj,j). This finding 
can be understood qualitatively from and ([53)1 . where one sees that r — ► /ry mainly 
depends on the large angle 623 while fi — > c-f involves the small quantities #13 and Am^. 
The difference in the scatter range of the predictions for r — > ^7 and \l — > thus reflects 
the different relative error of the quantities #235 #13 and Amj 2 and also the complete lack 
of knowledge on (p (see Tab. EI). Furthermore Figs. 3, 4 and 5 and Tab. El show that the 
experimental prospects favor the channel /i — > over r — > /X7 for testing small values of 
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Scenario 


Br(n -> e 7 ) = 1.2 • lO" 11 


Br(n -> e 7 ) = 10~ 14 


Br(r -> ^7) = 10~ 9 




M^/10 14 GeV 


Mr/10 


14 GeV 


M^/10 14 GeV 




present/hier. 


future/hier. 


future/deg. 


future/hier. 


A 


[4,20] 


[0.2,4] 


[0.9,3] 


- 


B 


[0.1,20] 


[0.006,0.4] 


[0.02,2] 


[1,2] 


C 


[0.6,20] 


[0.04,0.8] 


[0.2,2] 


[9,11] 


D 


[2,20] 


[0.07,2] 


[0.2,2] 


[15,19] 


E 


[0.3,20] 


[0.02,0.8] 


[0.1,2] 


[4,5] 


F 


[3,20] 


[0.2,2] 


[0.6,2] 


- 


G 


[0.2,20] 


[0.01,0.4] 


[0.1,2] 


[2,3] 


H 


[4,20] 


[0.3,4] 


[1,3] 


- 


I 


[0.04,5] 


[0.003,0.04] 


[0.02,1] 


[0.3,0.6] 


J 


[0.3,20] 


[0.02,0.8] 


[0.1,2] 


[3,4] 


K 


[0.5,20] 


[0.03,0.8] 


[0.2,2] 


[4,7] 


L 


[0.04,5] 


[0.003,0.04] 


[0.02,0.6] 


[0.2,0.5] 


M 


[0.6,20] 


[0.06,0.8] 


[0.2,2] 


[6,9] 



Table 3: Ranges of values for Mr for the given branching ratios in the case of hierarchical 
or degenerate neutrino spectra with present or future uncertainties of neutrino parameters. 
With '-' we denote that the sensitivity is too low. 



Mr. Larger values of Mr would be probed more accurately in r — > /i 7 . 
We also find that for fixed Mr the branching ratios for — > Zj 7 depend strongly on the 
particular mSUGRA scenario. The strongest bounds on Mr are obtained in scenario L 
due to very large tan f3 and small sparticle masses, whereas scenario H with large gaugino 
masses yields the weakest bounds. 

In summary, we find for hierarchical neutrino spectra that a future measurement of Br(r — > 
/i 7 ) 10~ 9 would typically determine Mr up to a factor of 2 given the uncertainties in the 
neutrino parameters. On the other hand, a measurement of Br(fi — > e 7 ) ~ 10~ 14 would 
determine the right-handed scale only up to a factor of 10-100, even if the SUSY parameters 
would be known. Finally, assuming an exactly massless lightest neutrino (as in previous 
works), the upper bounds on Mr improve by a factor of up to 10. For degenerate neutrinos 
and Mr < 10 14 GeV, Br(k — > ^7) is suppressed by roughly two orders of magnitude as com- 
pared to the case of hierarchical neutrino spectra, but exhibit a similar dependence on Mr. 
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This can be seen by comparing Figs. 3 and 4. However, for Mr > 10 14 GeV, the neutrino 
Yukawa couplings increase more strongly for degenerate than for hierarchical neutrinos, as 
illustrated in Fig. 2. Hence for sufficiently large M R) the branching ratios for hierarchical 
and degenerate neutrinos become comparable. This behaviour is particularly pronounced 
for Br{r — > /ry) as indicated in Fig. 5, because of the enhanced loop contribution from the 
lightest stau. 

6 Conclusions 

Future experiments searching for lepton-flavor violating rare processes can test the Majorana 
mass scale Mr of right-handed neutrinos in the see-saw mechanism. We have systematically 
and comprehensively studied the sensitivity of Br{li — > ^7) on Mr in mSUGRA benchmark 
scenarios designed for future collider studies taking into account the uncertainties of present 
and future neutrino measurements. We have assumed degenerate Majorana masses for 
the right-handed neutrinos and a normal neutrino mass hierarchy, and have considered 
hierarchical and degenerate neutrino spectra. 

For hierarchical neutrinos the measurement of Br\\i — * cy) ~ 10~ 14 would probe Mr in 
the range 5 ■ 10 12 GeV to 5 • 10 14 GeV, depending on the mSUGRA scenario. On the 
other hand, a future measurement of Br{r — * /I'j) at a level of 1CT 9 will determine M R in 
the range larger than 5 • 10 13 GeV with an accuracy of a factor of 2 for a given scenario. 
In the case of degenerate neutrino masses the upper bound on Mr which can be derived 
from Br(fi — > &y) < 10~ 14 is (1 — 3) • 10 14 GeV, independently of the mSUGRA scenario. 
Unification of the top Yukawa coupling and the Yukawa coupling of the heaviest neutrino 
at M x would fix Mr to M R pa 4 • 10 14 GeV and Mr pa 7 ■ 10 13 GeV for hierarchical and 
degenerate neutrinos, respectively. This proposition can thus be tested in the future. 
Planned measurements of ji — > 3e are not expected to improve the bounds on Mr. On 
the other hand, a future measurement of R {^i~Ti — > e~Ti) pa 10~ 18 is found to be more 
sensitive to Mr by a factor of about 2 than Br (// — > cy) ~ 10~ 15 . 

The correlation between the SUSY contribution 5a M to — 2 and Br(lj — > ^7) can be 
used to reduce the mSUGRA scenario dependence of the above tests. Comparing the ratio 
Br(r — * /i7)/(5a At ) 2 in Fig. 6 with the branching ratio for r — > 7x7 shown in Fig. 5, one can 
see that for fixed Mr and different mSUGRA scenarios the ratio varies by two orders of 
magnitude less than the value of Br(r fxj) itself. 
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Figure 2: Largest Yukawa coupling \Y U3 \ normalized to the top Yukawa coupling Y t at Mx 
for hierarchical and degenerate neutrino spectra (tan/3 = 30). The shaded area is excluded 

\Y I 2 

by the constraint " 3 < 0.3. 
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Figure 3: Branching ratio of \i — > &y for hierarchical neutrinos and uncertainties of future 
neutrino experiments in the mSUGRA scenarios leading to the largest (L, upper) and the 
smallest (H, lower) LFV rates. 
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Figure 4: Branching ratio of fi — > cy for degenerate neutrinos and uncertainties of future 
neutrino experiments in the mSUGRA scenarios leading to the largest (L, upper) and the 
smallest (H, lower) LFV rates. 
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Figure 5: Branching ratio of r — > /ry for hierarchical (upper) and degenerate (lower) neu- 
trino masses in the mSUGRA scenarios leading to the largest (L, triangles) and the smallest 
(H, circles) LFV rates. Open and filled symbols refer to neutrino measurements with present 
and future uncertainties, respectively. 
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Figure 6: Ratio Br[r — > fi^)/(5a^) 2 for hierarchical neutrinos and uncertainties of future 
neutrino experiments. Shown are the expectations for the mSUGRA scenarios F (upper) 
and C (lower) which embrace the predictions for all other benchmark scenarios of Tab. ^ 
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